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PREFACE 

The  work  described  in  this  report  was  authorized 
under  Project  No.  1L161101A91A,  Independent  Laboratory  In-House 
Research  (ILIR).  This  work  was  started  in  January  1981  and 
completed  in  December  1982. 

The  use  of  trade  names  or  manufacturers'  names  in 
this  report  does  not  constitute  an  official  endorsement  of  any 
commercial  products.  This  report  may  not  be  cited  for  purposes 
of  advertisement. 

Reproduction  of  this  document  in  whole  or  in  part  is 
prohibited  except  with  permission  of  the  Commander,  U.S.  Army 
Chemical  Research,  Development  and  Engineering  Center,  ATTN: 
SMCCR-SPS-T,  Aberdeen  Proving  Ground,  Maryland  21010-5423. 
However,  the  Defense  Technical  Information  Center  and  the 
National  Technical  Information  Service  are  authorized  to  reproduce 
the  document  for  U.S.  Government  purposes. 

This  report  has  been  approved  for  release  to  the 

public. 
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OPTICAL  CHARACTERIZATION  OF  TROPOSPHERIC  AEROSOLS 


1.  INTRODUCTION  AND  BACKGROUND 

With  the  increasing  availability  over  the  past  two 
decades  of  high-speed  electronic  computers,  complex  calculations 
of  optical  extinction  and  scattering  of  electromagnetic  radiation 
by  atmospheric  aerosols  using  the  Mie  theory!  have  become 
commonplace.  The  technology  has  advanced  so  swiftly  that 
specialists  measuring  cloud  extinction  of  aerosols  often  are 
unaware  of  advances  being  made  simultaneously  by  other  specialists 
measuring  angular  scattering  patterns  of  aerosols  with  nephelom- 
eters.  Such  measurements  and  calculations  can  be  done  in  numerous 
wavelength  regions,  but  the  technology  which  concentrated  first 
upon  the  visible  wavelengths  now  is  expanding  rapidly,  especially 
into  the  regions  of  relatively  high  atmospheric  transparency 
called  'windows'  in  the  infrared.  This  technological  expansion 
can  be  measured  by  a  sampling  of  papers  appearing  in  journals 
during  the  past  year,  discussing  aerosol  measurements  or 
calculations  of  extinction , 2 . 3 ,4 , 5 , 6 , 7 , 8 ,9  angular  light 
scattering  or  nephelometry , 5 » 10 ,11 ,12 , 1 3 ,14 , 15  iidar  16,17,18,19 
and  related  aerosol  technology  including  thermal  emissivity  20  an(j 
calibration  of  optical  particle  counters. 21 

Several  trends  have  been  noted  in  all  of  this  work. 
First,  the  utility  and  information  content  of  simple  extinction 
measurements  has  tended  to  be  overlooked  in  the  rush  to  exploit 
scattering  techniques.  Second,  the  polarization  techniques  do  not 
seem  to  include  applicable  techniques,  some  new,  that  can  be 
applied  both  to  extinction  and  to  polarization  measurements. 

Third,  if  taken  in  combination  at  this  point  in  time,  the  overall 
technology  might  be  further  advanced  than  is  realized  by  workers 
in  the  individual  disciplines.  For  example,  remote  characteri¬ 
zation  or  quantitative  analysis  of  tropospheric  aerosols  seems  a 
distinct  possibility,  both  by  conventional  transmission 
(extinction)  spectroscopy  and  by  angular  scattering  or  back- 
scattering  techniques.  This  report  examines  these  possibilities 
by  first  reviewing  the  technology  of  aerosol  extinction 
spectroscopy. 

The  Beer-Lambert  equation  for  atmospheric  aerosols  can 
be  written: 

(T\)  =  a\  CL  (1) 


where  for  a  given  wavelength  \(nm) ,  T\  is  the  fractional  trans¬ 
mittance,  is  the  optical  mass  extinction  coefficient  of  the 
aerosol  (m2g-l),  C  is  the  aerosol  mass  concentration  (gm~3)  and 


9 


WTWWfWW 


w 


W1K’WW',AM . V  7V.“  *7  '.V' V  «»  TO  »>U  VVW  TOT  ■ 


I 


L  is  the  optical  path  length  (m).  Some  form  of  Equation  1  has 
been  used  by  spectroscopists  virtually  since  the  birth  of  spec¬ 
troscopy  and  the  equation  has  become  so  well  known  that  the  power 
of  simple  extinction  measurements  of  aerosols  often  is  not  fully 
appreciated.  For  example,  using  suitable  multi-wavelength 
techniques,  extinction  measurements  can  contain  sufficient 
information  to  allow  direct  determination  of  the  mean  size  of  a 
spherical  droplet  aerosol  distribution  or  even  remotely  to 
characterize  chemical  reactions  in  cloud  droplets  at  isosbestic 
wavelengths,  without  the  use  of  polarized  radiation  or  other  more 
sophi sti cated  techniques .22,23,24,25 

It  is  not  surprising,  however,  that  polarization 
techniques  have  received  heavy  emphasis  since  their  exploitation 
in  the  early  1940's  to  characterize  military  fog-oil  smoke  mean 
droplet  diameter  based,  in  part,  in  earlier  work  at  the  U.S. 

Bureau  of  Standards  and  by  Stratton  and  Houghton  using  the  Mie 
theory. 26, i  The  polarization  ration  then  was  measured  only  for 
transparent  spheres  for  which,  if  the  complex  index  of  refraction 
m^  is  given  by  (n  -  ik)x  where  is  the  real  index  and 
k^  is  the  imaginary  index  (closely  related  to  the  absorption 
coefficient),  then  k^  =  0.  Polarization  measurements  can  be  used 
to  determine  particle  size  distribution  in  addition  to  mean  size 
of  aerosol  droplet  distributions.  This  was  recognized  and  the 
technique  was  developed  in  the  1960 ' s . 27 , 28 , 29  gy  1959 y  the 
techniques  of  particle  size  distribution  measurements  had  been 
extended,  as  discussed  in  Kerker's  classical  book  on  light 
scattering.  This  book  gives  an  excellent  review  of  the  technology 
as  it  existed  then. 20 

In  recent  years,  there  is  an  increased  awareness  that 
aerosol  particle  absorption,  e.g.,  in  the  infrared  where  k  ^  o, 
can  be  interpreted  in  new  ways  utilizing  not  only  extinction 
measurements  or  polarization  measurements,  but  the  two 
together. 22, 23, 24, 30  The  author  concurs  that  this  is  a  promising 
field  for  new  research  and  this  Laboratory  is  developing  an 
expanded  capability  for  scattering  and  polarization  measurements 
to  complement  ongoing  research  in  extinction  spectroscopy  of 
aerosols . 


2.  AEROSOL  EXTINCTION  MEASUREMENTS 

It  is  known  that  the  variation  of  k^  from  values  near 
zero  in  the  visible  wavelengths  to  significant  ones  in  the 
infrared  wavelengths  provides  a  method  for  the  determination  of 
aerosol  mean  droplet  size  and  mass  concentration  from  simple 
extinction  measurements  if  two  or  more  wavelengths  are  selected 
properly  for  observation.  For  example,  at  the  He:Ne  wavelength 
X  =  0.63  4m,k^0  and  typical  water  fog  droplets  have  diameters  D4 
(in  4m)  such  that  D4  » ^ or ,  to  use  Kerker's  criterion,  the  size 
parameter  ttD^  /  X  >  2.0  so  that  the  conditions  for  geometric 
scattering  exist.  This  leads  to  a  constancy  of  the  product 
3. 0-3. 2,  where  is  defined  in  Equation  1.  At  the  same  time, 
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the  product  a ^  •  Dp  decidedly  is  not  constant  with  Dp  at  wave¬ 
lengths  wh«..re  k^o  in  the  infrared.  This  is  illustrated  in  Figure 
1,  calculated  from  the  Mie  theory.  In  Figure  1,  curves  are  shown 
for  the  infrared  wavelengths  of  8.5,  10.5  and  12.57  (or  4?r)  pm,  as 
well  as  for  a  composite  of  wavelengths  in  the  8.5  -  12.57  pm 
region  where  extinction  can  be  measured  by  a  simple  broadband 
tratismissometer  operating  at  8-13  pm  in  the  atmospheric  window 
region.  While  operation  at  a  specific  wavelength  such  as  the  10.6 
pm  CO2  laser  line  for  comparison  to  a  visible  wavelength  such  as 
the  He:Ne  line  at  X=  0.63  pm  has  advantages,  a  broadband  infrared 
reference  band  also  can  be  used.  In  fact,  related  findings,  using 
a  filtered  light  source  in  the  visible  at  X=  0.515  pm  and  a  9-12 
pm  broadband  transmissometer  to  study  the  formation  and 


were  reported  at  least  as  early  as 


When  the  techniques  of  Carlon  et  al.,22,24  are  extended 
to  use  the  Mie  theory  to  calculate  functions  of  the  transmittances 
at  =  0.63  p m  and  the  broadband  or  composite  wavelength  band 
8.5-12.57  pm  versus  water  droplet  diameter  Du,  curves  like  those 
shown  in  Figure  2  are  obtained.  It  is  seen  that  the  calculated 
curves  are  widely  separated.  Thus,  a  sensitive  method  is  provided 
for  the  determination  of  mean  droplet  diameter.  Experimental 
points  for  a  steam-generated  water  fog  are  shown  in  Figure  2. 

These  give  excellent  agreement  with  calculated  curves,  for  example 
Dp.  =  8  pm.  This  result  recently  was  cited  by  Weinman  et  al.,^2  as 
approximating  Deirmend  jian'  s^3  mcfdel  C-l  (or  vice  versa),  where  Dp 
is  taken  as  the  mass  median  diameter.  These  results  seem 
credible.  What  is  uncertain  is  the  kind  of  distribution  to  which 
this  diameter  is  most  applicable.  Kerker^O  has  given  an 
informative  discussion  of  particle  size  distributions,  and  of 
transmission  (extinction)  and  turbidity  measurements.  In  this 
report,  turbidity  can  be  taken  as  the  product  ay C  in  Equation  1, 
having  the  units  m~l . 


Justification  for  the  use  of  composite  of  wavelengths 
(Figure  1  and  abscissa,  Figure  2)  can  be  found  by  examining 
calculated  curves  and  steam  generated  fog  data  points  for  the 
individual  wavelengths  X=  8.5,  10.5  and  12.57  (or  bir)  pm  and 
plotting  curves  like  Figure  2  as  is  done  in  Figure  3.  Viewing  the 
plots  from  top  to  bottom  in  Figure  3,  it  can  be  seen  that  the  fog 
data  gave  best  agreement  with  calculation  for  a  mass  median 
diameter  (MMD)  of  6  pm  at  X  =  8.5pm,  MMD  =  10  pm  at  X  =  10.5  pm 
and  MMD  =  6  pm  at  X  =  4  pm.  In  his  discussions  of  particle  size 
distribution  functions  and  transmission  measurements,  Kerker^O 
cautions  that  mean  particle  diameters  obtained  using  methods  like 
those  discussed  in  the  present  paper  may  be  adequate  for  most 
purposes  if  the  distribution  is  sufficiently  narrow,  but  strictly 
are  valid  only  for  monodisperse  systems.  By  using  a  composite  of 
wavelengths  to  obtain  a  mean  diameter,  one  is  averaging  both  over 
wavelength  intervals  and  over  particle  size  distributions.  The 
results  indicate  that  the  technique  works  very  well,  especially 
under  dynamic  droplet  growth  or  evaporation  conditions  where  what 


Figure  1.  Values  of  aX  D/a  vs.  D/a  Calculated  from  the 

Mie  Program  for  Water  Using  the  Data  of  Hale 
and  Querry.38  The  Average  Value  of  the 
'Constant'  Curve  Is  Shown;  Actually  the 
Function  Is  a  Damped  Oscillation. 
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Figure  2. 


Comparison  of  Computed  Curves  Relating 
Visible  to  Infrared  Optical  Transmittance 
of  Water  Fog  to  Experimental  Data  Points  for 
Cooling,  Steam-Generated  Water  Fog  Clouds; 
Note  that  Data  Points  Cluster  Along  the  Curve 
Dp  =  8  p.m  in  Good  Agreement  with  Gravimetric 
Data  from  the  Same  Experiment. 
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At  wavelengths  called  'isosbestic  points', 25  the  complex 
indices  (n  -  ik)\  of  some  liquids  are  such  that  remains  con¬ 
stant  even  though  the  liquid  solutions  comprising  the  droplet 
aerosol  vary  widely  in  their  chemical  composition,  e.g.,  for 
droplets  of  orthophosphoric  and  related  acids  that  might  be  formed 
by  burning  phosphorus  in  air  of  varying  relative  humidity  such 
that  the  solute  (acid)  concentrations  vary  from  near  zero  to  near 
8 5%  by  weight.  In  Figure  5,  X  =  11.4  ^.m  is  an  isosbestic  point 
calculated  from  refractive  index  measurements  of  pure  ortho¬ 
phosphoric  acid  (H3PO4)  in  water. 36 » 37  Real  atmospheric  smokes, 
obtained  by  burning  phosphorus,  may  contain  mixtures  of  acids 
depending  on  the  rates  of  combustion  in  the  presence  of  varying 
humidities,  and  so  may  not  give  experimental  spectra  like  those 
for  pure  H3PO4  in  Figure  5.  where  a  11,4  remains  essentially 
constant  at  0.13-0.15  m2g-i,  while  acid  concentration  varies  from 
near  zero  to  85%  by  weight.  This  will  be  investigated.  Because 
of  the  constancy  of  a  11.4  compared  to  the  wide  range  of  values  of, 
say  a  9.7  (see  Figure  5)  with  solute  concentration,  remote 
characterization  of  the  aerosol  becomes  a  possibility  if  something 
is  known  about  the  droplet  constituents.  It  can  be  shown  that: 

**  4nk^f  (m^)/Ap  (  2 ) 


where  p  (g  cm-3)  is  the  droplet  solution  (acid)  density  i.e.,  mass 
density,  and: 


Figure  5.  Computed  Values  of  a\  for  Liquid  Droplet 
Aerosols  Comprising  Several  Concentrations 
of  Orthophosphoric  Acid  (H3PO4)  in  Water, 
from  Optical  Constants  of  Querry  and  Tyler.  37 
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where  m\  =  (n  -  ik)\  and  the  function  of  Equation  3  is  plotted  in 
Figure  6  for  values  of  rj ^  marked  on  the  curves  and  ranging  from 
0.1  to  4.0,  and  for  values  for  k\  shown  on  the  abscissa.  It  is 
emphasized  that  the  above  discussion  applies  to  Rayleigh- 
scattering  smokes  where  D^.  <<^  as  clearly  is  the  case  for 
orthophosphoric  or  other  acid  aerosols  (smokes)  produced  by 
burning  phosphorus  in  moist  air. 

Often  it  is  useful  to  plot  extinction  coefficients  or 
fractional  transmittances  versus  Dp  rather  than  versus  the  size 
parameter  -tr  D ^.1  \  t  for  specific  wavelengths.  Two  such  calculated 
plots  are  shown  for  water  droplets  in  Figure  7  and  8.38  jn  Figure 
7,  the  curves  are  calculated  for  a  product  C  x  L  =  5.0  (Equation 
1)  and  for  wavelengths,  reading  from  top  to  bottom  on  the  left- 
hand  ends  of  the  curves,  of  X  =  0.63,  8.5,  10.5  and  4  rr  pm  (the 
same  wavelengths  considered  in  Figures  1-3).  In  Figure  8,  the 
curves  are  calculated  for  X  =  10  pm  to  show  the  functions  vs.  D^. 
of  total  extinction,  and  of  the  scattered  and  absorbed  components 
that  are  summed  to  give  total  extinction,  i.e., 

aT=as  =  +aA  (4) 

where  the  subscripts  indicate  total  (T) ,  scattering  (S)  and 
absorptive  (A)  extinction. 

Sassen23  recognizes  that  aerosol  cloud  extinction 
measured  in  conjunction  with  angular  scattering  measurements  could 
lead  to  remote  sensing  of  cloud  composition  by  using  two  wave¬ 
lengths  like  X  =0.63  and  10.6  pm.^  But  the  present  discussion 
indicates  that  such  remote  characterization  could  solely  result 
from  cloud  extinction  measurements  at  two  closely-spaced  infrared 
wavelengths,  one  of  them  an  isosbestic  point.  Additional  data 
might  be  obtained  by  combining  extinction  and  angular  scattering 
measurements  at  an  isosbestic  wavelength.  Optical  parameters  can 
carry  only  so  much  information,  but  the  idea  is  provocative  enough 
to  warrant  an  investigation. 

Other  investigations  for  remote  characterization  of 
tropospheric  aerosol  clouds  are  possible,  such  as  those  based  on 
effects  peculiar  to  Christiansen  wavelengths . 39  a  complete 
disctjssion  is  given  by  Carlon  where  the  latter  reference  discusses 
aerosols  ranging  in  size  from  molecular  clusters  to  water  cloud 
droplets .  2  >  39 

3.  POLARIZATION  AND  ANGULAR  SCATTERING  MEASUREMENTS 

Kerker30  has  given  an  excellent  review  of  techniques  to 
determine  particle  size  distributions  using  these  measurements  and 
I  will  not  elaborate  except  to  raise  a  few  pertinent  points. 
Sassen23  has  verified  that  angular  scattering  patterns  calculated 
from  the  Mie  theory  can  be  verified  rather  closely  by  experiment 
not  only  for  spherical  particles  but  for  irregularly-shaped  ones 
as  well.  Thus  polarization  or  angular  scattering  techniques 
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Figure  7.  Fractional  Transmittance  for  a  C  x  L  Product 
=  5.0  (Equation  1)  of  Water  Fog  vs.  Droplet 
Diameter  for  Wavelengths,  Reading  from  Top 
to  Bottom  of  the  Left-Hand  Ends  of  the  Curves, 
of  x  =  0.63,  8.5,  10.5,  and  kir  *im,  which  Are 
the  Wavelengths  also  Considered  in  Figures 
1-3. 


Figure  8.  ax  vs.  Calculated  from  Mie  Program 
Using  Data  of  Hale  and  Querry38  at 
X  =  10  #i.m  for  Water  Droplets;  the 
Subscripts  T,  S,  and  A  Refer  to  Equation 
(4)  and  Definitions. 
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combined  with  some  of  the  recent  advances  in  extinction  measure¬ 
ment  technology  discussed  here  should  provide  new  combinations  for 
the  remote  characterization  of  tropospheric  aerosols. 

The  polarization  method  for  evaluating  the  distribution 
of  sphere  sizes,  utilizes  polarization  of  scattered  radiation 
obtained  from  a  monochromatic  light  source  at  various  angles  of 
observation . 30  The  intensity  of  the  component  of  the  scattered 
light,  whose  electric  vector  vibrates  perpendicular  and  parallel 
to  the  plane  of  observation,  1^(0)  and  12(d) ,  is  measured  at  a 
number  of  angles,  where  6  is  the  scattering  angle  measured  from  the 
direction  of  the  incident  light  beam.  A  polarization  ratio,  P(0), 
is  defined  such  that: 

P(0)  =  I2  (0)  /  Tx  (9)  (5) 

For  example,  P ( 0 )  can  be  calculated  as  a  function  of 
for  any  assumed  particle  distribution  such  as  the  zeroth  order 
logarthmic  distribution  (ZOLD)  and  for  any  standard  deviation, 

0o,  ranging  from  0o  =  0  (i.e.,  a  monodisperse  aerosol),  to  values 
of  0.3  or  more  where  the  curves  flatten  and  the  method  approaches 
its  limit  of  usefulness.  Figure  9  shows  the  family  of  curves 
calculated  for  the  He:Cd  wavelength  X  =  0.4416  pim,  =  0.3  ^m, 
i)X=  1.484  and  k  \  =  0  so  that  m^  =  (1.484  -  i •  o )  \  and  the  density 

is  p  =  0.98  g  cm  "3.  The  curves  are  for  dioctyl  phthalate  (DOP),  a 
liquid  commonly  used  for  test  aerosols. 2  Agreement  between 
calculation  and  measurement  using  the  polarization  method  is  good 
to  excellent,  often  within  a  few  percent. 

The  effect  of  the  wavelength  of  the  monochromatic 
light  upon  the  function  shown  in  Figure  9  is  of  special  interest. 
This  has  been  investigated^ •  > 28  for  visible  and  near  infrared 
wavelength,  including  cases  such  as  vanadium  pentoxide  (V2O5)  for 
which  strong  absorption  occurs  at  shorter  visible  wavelengths, 
i.e.,k\*0.Z9  Special  techniques  using  extinction  measurements  in 
the  infrared  where  k^*0,  suggest  that  combinations  of  visible  and 
infrared  wavelengths  can  optimize  desired  parameters  using 
polarization  techniques.  Since  the  size  parameter  tt  D/a  /  X 
depends  on  wavelength,  a  change  in  the  wavelength  of  the  incident 
light  will  cause  a  shift  in  the  position  of  the  curve  peaks  in 
Figure  9,  larger  values  of  X  shift  the  peaks  to  the  right. 
Multi-wavelength  observations  have  been  used  to  show  that  particle 
size  distributions  obtained  at  different  wavelengths  are 
consistent,  but  when  the  wavelength  spacing  is  as  great  as  X  =  0.63 
fi.m  and  10.5  /i.m,  additional  work  is  needed  to  confirm  this  and  also 
to  determine  whether  the  useful  range  of  0o  can  be  extended  beyond 
0.3  in  measurements  using  the  polarization  method. 30 

The  scattering  ratio  method  for  the  determination  of 
particle  size  distribution  is  a  variation  where  the  angle  of 
observation  is  held  constant  and  the  polarization  ratio  is 


measured  or  calculated  over  a  spectrum  or  over  a  number  of 
discrete  wavelengths  .40, 41 ,42  The  method  was  perfected  by  using 
monodlsperse ,  polystyrene,  latex  spheres  and  known  combinations  of 
single  sphere  sizes  to  give  mixtures  of  known  size  distributions 
checked  by  electron  microscopy.  The  scattering  ratio  method  gives 
an  alternative  way  of  studying  the  variation  of  0o  at  some  D fj.  with 
X  at  some  constant  6  .  The  same  thing  can  be  done  by  using  the 
polarization  method  directly  and  plotting,  e.g.,  Figure  9  with  X 
as  the  abscissa  at  constant  6. 

Turbidimetric  particle  size  distribution  methods 
invariably  are  similar  to  those  of  Wallach  et  al.40  Very  recently 
Nelson43  proposed  a  method  by  which  calculated  scattering  cross 
sections  from  the  Mie  theory  are  compared  to  experimental 
extinction  data  using  monodi spersed ,  polystyrene,  latex  spheres  in 
distilled  water  as  a  calibration  means.  A  computer  program  is 
used  to  null  the  error  versus  wavelength,  forcing  a  unique  solution 
for  the  size  distribution  which  is  taken  as  bimodal  to  account  for 
coagulation.  The  method  also  is  said  to  yield  particle  volume 
concentration  and  the  scattering  cross  section  function,  (Q)\. 
Extrapolation  of  the  technique  to  tropospheric  aerosols  remains  to 
be  demonstrated. 

4.  SUMMARY  AND  CONCLUSIONS 

This  report  has  concentrated  on  advances  that  have  been 
made  in  the  remote  characterization  of  tropospheric  aerosols  by 
extinction  and  turbidimetric  methods  over  the  decade  that  has 
elapsed  since  the  publication  of  Kerker's  classic  text. 30  The 
technology  has  been  updated,  including  results  not  previously 
published,  with  improved  techniques  in  polarization  and  angular 
scattering  methods.  Perhaps  most  importantly,  the  author  has 
tried  to  support  the  argument  that  if  all  methods  presently 
available  in  this  technology  are  taken  into  account  at  this  time, 
our  ability  to  remotely  characterize  or  quantitatively  analyze 
atmospheric  aerosols  is  further  advanced  than  workers  in  the 
individual  light  measurement  disciplines  realize. 

It  is  concluded  that  by  combinations  of  well  established 
and  recently  developed  light  scattering  techniques,  investigators 
can  undertake  promising  new  areas  of  research.  Some  examples 
include  combinations  of  extinction  and  polarization  measurements 
with  computer  programs  to  study  isosbestic  or  Christiansen 
wavelengths,  and  optimization  of  desired  parameters  as  functions 
of  wavelength  or  particle  size  distribution. 
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